Time-resolved measurement of Landau— Zener tunneling in different bases 
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A comprehensive study of the tunneling dynamics of a Bose-Einstein condensate in a tilted 
periodic potential is presented. We report numerical and experimental results on time-resolved 
measurements of the Landau-Zener tunneling of ultracold atoms introduced by the tilt, which 
experimentally is realized by accelerating the lattice. The use of different protocols enables us to 
access the tunneling probability, numerically as well as experimentally, in two different bases, namely, 
the adiabatic basis and the diabatic basis. The adiabatic basis corresponds to the eigenstates of 
the lattice, and the diabatic one to the free-particle momentum eigenstates. Our numerical and 
experimental results are compared with existing two-state Landau-Zener models. 
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I. INTRODUCTION 

Quantum transport is an essential topic in solid state 
physics and electronic applications. Bloch oscillations, 
Landau-Zener (LZ) tunneling, and Wannier-Stark lad- 
ders [l|-|9|, are fundamental quantum effects occurring 
in a system of electrons moving in a periodic potential 
and driven by an electric field. Due to complications 
such as impurities, lattice vibrations, and multiparticle 
interactions, clean observations of these effects have been 
difficult [lOj . In recent years, ultracold atoms and Bose- 
Einstein condensates in optical lattices have been increas- 
ingly used to simulate solid state systems and the above 
mentioned phenomena [aTS IflHla ] . 

Optical lattices are easy to realize in the laboratory, 
and their parameters can be perfectly controlled both 
statically and dynamically, which makes them attractive 
as model systems for crystal lattices. The LZ model for 
transitions 0, y] between two energy states at an avoided 
level crossing is one of the few exactly solvable examples 
of time-dependent quantum mechanics. LZ transitions 
have been investigated for Rydberg atoms [17|, molecu- 
lar nanomagnets [13, UM i field-driven superlattices [2(| , 
current-driven Josephson junctions [2l| , Cooper-pair box 
qubits [2211 . and using light waves in coupled waveg- 
uides [2a - l25j . While the asymptotic tunneling probabil- 
ity can be calculated accurately [26[ and has an intuitive 
interpretation as a statistical mean value of experimental 
outcomes, the concept of tunneling time and its compu- 
tation are still the subject of debate even for simple sys- 
tems [2l|, l27H3l| . The tunneling time is the time required 
for a state to evolve into an orthogonal state. 

In this paper, we present numerical as well as experi- 
mental results on the Wannier-Stark system. This sys- 
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tern is realized with ultracold atoms, forming a Bose- 
Einstein condensate, in an optical lattice subjected to a 
static tilting force [8|]. The tilt is experimentally imple- 
mented by accelerating the optical lattice (q.M. ll3 . [lal32| - 
l36j |. We explore the LZ tunneling between the Bloch 
bands of a Bose-Einstein condensate in such an accel- 
erated lattice. The lattice depth controls the tunneling 
barrier, while its acceleration controls the time depen- 
dence of the Hamiltonian. At large accelerations LZ tun- 
neling leads to significant interband transitions for the 
condensate [15|, [26[ . This tunneling process is detected 
by measuring the atomic momentum distributions. 

Following our previous work, in which we presented 
time-resolved observations of LZ tunneling |36| . in the 
present article we report more detailed investigations. 
We measure the time dependence of the tunneling proba- 
bility by performing a projective quantum measurement 
on the eigenstates in a given basis of the Hamiltonian de- 
scribing the Bose-Einstein condensate within the optical 
lattice. Our measurements resolve the steplike time de- 
pendence of the occupation probability. Using different 
numerical as well as experimental protocols, we are able 
to perform our calculations and experiments both in the 
adiabatic basis of the lattice eigenstates and in the di- 
abatic basis of the free-particle momentum eigenstates. 
We present theoretical and experimental results which 
clearly show that the time dependence of the transition 
probability exhibits a stcplike structure with a finite tran- 
sition time and oscillations with a finite damping time, 
all of them depending on the choice of the measurement 
basis. To our knowledge, such time-resolved measure- 
ments in different bases have not been reported for other 
systems before. 

The paper is organized as follows. Section II col- 
lects the necessary theoretical background to describe 
the probability and transition time for the LZ transition 
tunneling. The limits one faces in applying this theory to 
the Wannier-Stark problem, and the essential theoretical 



and numerical tools to describe our time-resolved mea- 
surements are reported in Section III. Section IV presents 
numerical and experimental data. We discuss and sum- 
marize our results in Section V. 



II. SURVIVAL PROBABILITY AND 
TRANSITION TIME 

A. LZ theory in a nutshell 

Quantum mechanical systems having two discrete en- 
ergy levels arc omnipresent in nature. For crossing lev- 
els, there is a possibility of a transition if the degeneracy 
is lifted by a coupling and the system is forced across 
the resulting avoided crossing by varying the parameter 
that determines the level separation. This phenomenon 
is known as a LZ tunneling. LZ theory, developed in 
the early 1930's in the context of atomic scattering pro- 
cesses and spin dynamics in time-dependent fields [l|- 
|j], demonstrated that transitions are possible between 
two approaching levels as a control parameter is swept 
across the point of minimum energy separation. The 
phase accumulated between the incoming and outgoing 
passages varies with, e.g., the collision energy, giving rise 
to Stuckelberg oscillations in the populations |3|. 

In its basic form the LZ problem can be described by a 
simple two-state model and allows for a simple expression 
for the transition probability. The LZ Hamiltonian for a 
single crossing taking place at time t = can be written 
as the following 2 by 2 matrix 



(at AE/2 \ 
Hlz - { AE/2 -at ) 



The diabatic basis is the eigenbasis of the bare states 
of Eq. (TTJ) when there is no off-diagonal coupling. The 
adiabatic basis, on the other hand, is the basis of a 
system with a finite AE/2 coupling between the two 
states. The Hamiltonian has two adiabatic energy lev- 



els E± = ±\\J{2cd) + AE 2 . 

Assuming that the system is initially, at t — > — oo, in 
the ground energy level E- and if the sweeping rate is 
small enough, it will be exponentially likely that the sys- 
tem remains in its adiabatic ground state E— at t — > +oo. 
The limiting value of the adiabatic LZ survival probabil- 
ities (for t going from — oo to +oo) is |26j . 



P a (oo) = 1 - cxp 



7 



(2) 



(1) 



where we introduce a dimcnsionlcss parameter, the so 
called adiabaticity parameter 7 = Aha/ AE 2 . This sur- 
vival probability is valid for both 25_ and E+ initial 
states, and the same equation is valid for the diabatic 
case. A small adiabaticity parameter corresponds to a 
small velocity of the state displacement along the en- 
ergy scale compared to AE 2 , such that the system fol- 
lows the adiabatic trajectory of Fig. Q] Thus, there is 
a large coupling between the diabatic states and at the 
avoided crossing at t = an almost complete transition 
from the initial diabatic state to the final diabatic state 
takes place. On the other hand, for a large value of the 
adiabaticity parameter 7, the coupling between the two 
states is small and consequently the system remains in its 
initial state following the diabatic trajectories of Fig. [1] 



B. Jump times 




2at/AE 

FIG. 1: Energy levels as a function of time. The dashed lines 
show the so-called diabatic levels, i.e., the energies of states in 
the absence of the interaction. The solid lines demonstrate the 
so-called adiabatic levels, i.e., the eigenstates of the system 
corresponding to the instantaneous Hamiltonian. 



The off-diagonal term, AE/2, is the coupling between 
the two states and a is the rate of change of the en- 
ergy levels in time. The dynamics of the system can 
be measured in different bases, diabatic and adiabatic. 



A careful study of the transition from an initial state 
to a final state can reveal the time required to complete 
the transition. Moreover, in the case of multiple level 
crossings, as in our experimental realization of ref. [361 ] . 
it is necessary to know whether a transition has been 
completed before the next avoided crossing. The LZ ap- 
proach may be applied when a transition between two 
coupled quantum states takes place in a small time inter- 
val around the avoided crossing and successive crossings 
are independent from each other. 

Analytical estimates for the LZ transition times have 
been derived in [29|, |3fJ using the two-state model of 
Eq. ([1]). In a given basis, e.g., adiabatic or diabatic, 
different transition times are obtained. Vitanov |29j cal- 
culated the time-dependent diabatic/adiabatic survival 
probability at finite times. The LZ transition times were 
derived in [30j using some exact and approximate results 
for the transition probability. Fig.[2]shows a typical time 
dependence of the adiabatic survival probability, similar 
to that predicted in [29|], that we measured for Bose- 
Einstein condensates in optical lattices for experimental 
parameters to be discussed in Section llVl Notice that in 
the Bose-Einstein condensate case the crossing occurs at 



the time t — Tb/2 defined below. The t 
value is given by Eq. ©. 
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FIG. 2: (Color online) Time-resolved measurements of the 
adiabatic survival probability for the Bose-Einstein conden- 
sate tunneling in an optical lattice at fixed dimensionless force 
F = 1.197 and different lattice depths: V/E Iec = 2.3 (filled 
squares), 1.8 (open circles), 1 (open squares) and 0.6 (filled 
circles); all parameters are introduced in detail in section IIIA. 
For the Bose-Einstein condensate evolution the crossing time 
is t — Tbioc1i/2, where the step of the survival probability is 
one half of the final value. The dashed and dot-dashed lines 
are the results of numerical simulations using the cut-off and 
adiabatic method, respectively (see section IIIC). The lattice 
depth for the numerical simulations was corrected by up to 
±15% with respect to the experimentally measured values to 
give the best possible agreement. 



The LZ jump time in a given basis can be defined as 
the time after which the transition probability reaches its 
asymptotic value. From this definition one can expect to 
observe a step-like structure, with a finite width, in the 
time-resolved tunneling probability, as in Fig. [2l Because 
the step is not very sharp, it is not straightforward to 
define the initial and final times for the transition. It is 
even less obvious how to define the jump time for both 
small and for large coupling. Some possible choices have 
been used by Lim and Berry |28( and Vitanov [29|, [3(| ■ 
The problem is even more complicated when the survival 
probability shows an oscillatory behavior on top of the 
step structure as seen in Fig. [21 which shows experimental 
and numerical results for a single LZ transition measured 
in the adiabatic basis (the numerical and experimental 
methods will be described in detail later in this paper). 
The oscillations give rise to other time scales in the sys- 
tem, namely an oscillation time and a damping time of 
the oscillations appearing in the transition probability 
after the crossing. Therefore, a measurement of the tun- 
neling time depends very much on how these times are 
defined and also which basis is considered. 



In [29| , the jump time in the diabatic/adiabatic bases 
is defined as 



jump d/ 

r d/a 



i(oo) 



^/a(0) ' 



(3) 



where -Pd/a is the transition probability between the two 
diabatic/adiabatic states, respectively. i^/ a (0) denotes 
the time derivative of the tunneling probability evalu- 
ated at the crossing point. The diabatic jump time 
T ramp ^ ^/2irh/a is almost constant for large values of 



the adiabaticity parameter 7 [2S 
decreases with increasing 7, T^ um ' 



Instead, for 7 <C 1 it 
~2$ 



2 x A( 7 a)- 1 [29|. In 
the adiabatic basis, when 7 is large the transition prob- 
ability resembles the one of the diabatic basis with an 
equal jump time. For a small adiabaticity parameter, 
because of the oscillations appearing on top of the tran- 
sition probability step structure, it is not straightforward 
to define the initial and the final time for the transition. 
Vitanov [29j defines the initial jump time as the time 
t < at which the transition probability is very small 
(i.e., Ph(t) = eP a (oo), where e is a proper small num- 
ber). The final time of the transition t > is defined 
as the time at which the non-oscillatory part of P a (r) 
is equal to (1 + e)P a (oo). Using these definitions, Vi- 
tanov derived that the transition time in the adiabatic 
basis depends exponentially on the adiabaticity parame- 
ter, tI u " 1 p w (4/e) 1/6 7~ 1/3 cxp(7r/(6 7 )) y/Ji/a, |tg|. 

In principle, the experimental and numerical methods 
presented in the following could be used for a quantitative 
study of the tunneling time (or jump time) as function 
of the parameters of the system. For the purposes of the 
present paper, however, we concentrate on a careful anal- 
ysis of the possibilities and limitations of our methods, 
and in particular on measurements of LZ tunneling in 
different bases. 



III. LZ IN AN OPTICAL LATTICE POTENTIAL 

A. Wannier— Stark problem and LZ limit 

We generalize the two-level LZ theory to study the 
temporal evolution of ultracold atoms loaded into a spa- 
tially periodic potential subjected to an additional static 
force in the presence of negligible atom-atom interac- 
tions, as in the experimental conditions |36|. The dy- 
namics of ultracold atoms in a tilted optical lattice can 
be described by the well-known Wannier-Stark Hamilto- 
nian |9| 



H 



2M dx 2 



V 



cos(2k L x) + F LZ x, 



(4) 



where M is the atomic mass, V is the depth of the optical 
lattice with the spatial period g?l = Al/2 determined by 
the laser wavelength Al, ^l = 27t/Al is the wave number 
of the laser light creating the periodic potential, and F^z 



is the Stark force. The characteristic energy scale of the 
system is the recoil energy which is defined as E ICC = 
n 2 h 2 /2Mdl. 



written in dimcnsionlcss units [26l . 137 



2dxl 



ff = -^TJ-2 + T7r COs(x ) 
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FIG. 3: (Color online) Band structure in the optical lattice 
potential and experimental protocols for measuring LZ dy- 
namics in the adiabatic and diabatic bases. After the initial 
loading into the lattice and acceleration for a time £lz (top), 
measurements of the instantaneous populations in the two 
states are performed (bottom) as explained in the text. In 
the top figures, the adiabatic (solid lines) and diabatic ener- 
gies (dashed lines) for an optical lattice of depth Vo — 1 are 
shown. 



The atomic motion produced by the force -FLz may be 
interpreted in the upper left energy diagram of Fig. [3] 
where for the case of Flz = the atomic energies E(q) 
for the n = 0,1,2 lower bands are represented versus 
the quasimomentum q within the Brillouin zone of width 
Pb = 2pio C = 2-Kfi/d^ [13, l32| ■ Under the action of a 
constant force -FLz , the quasimomentum of a condensate 
initially prepared at q = in the n — band scans 
the lower band in an oscillating motion periodically with 
the Bloch period T B = 2h(F LZ d L )~ 1 . At the edge of 
the Brillouin zone, where a level splitting AE increasing 
with V [26| takes place, tunneling of the condensate to 
the n = 1 energy band may occur. 

The Wannier-Stark Hamiltonian of Eq. @ can be 



where xq = 2iTx/di J , and energy and time are rescalcd 
by H Q = H/(8E lec ) and t = 8tE rcc /h respectively. 
Moreover, in dimcnsionlcss units, the lattice depth 
is given by Vo = V/E Tec and the force by F = 
Fhzdh/E Tec . The translational symmetry of the given 
Hamiltonian, broken by the static force, is recovered us- 
ing a gauge transformation. Substituting ip(xo,to) = 
cxp(— iFQtQXo/lQi:)%l)(xQ, to), the Schrddinger equation 
reads idtp/dto = H(t )ip, with H (to) the time-dependent 
Hamiltonian 



H(t ) 



P 



16tt 



Vo_ 
16 



cos(x ), 



(6) 



and the momentum operator p = —id/dxQ. In the fol- 
lowing, we analyze the Hamiltonian of Eq. (|6]) in the mo- 
mentum basis. In order to decompose the Hilbcrt space 
into independent subspaces, we use the Bloch decompo- 
sition and for that we identify the momentum cigenstatcs 
of the free particle (Vq = = Fq) for fixed quasimomen- 
tum q within the Brillouin zone, i.e., p = q + n, p and 
q being indices in the momentum and quasimomentum 
representations and n € Z. To calculate the time evolu- 
tion of any momentum eigenstate \p) = \q + n), we only 
need the Hamiltonian H q acting on the subspace with a 
given quasi-momentum index q, as there is no transition 
between states with different q 



1 



H -=2 



{q-lf Vo/16 
V16 (v) 2 Vb/16 
Vb/16 (q+1) 2 







\ 



V o 



(7) 



/ 



where q = q — Fo£o/167r. 

The full dynamics of the Wannier-Stark system can be 
locally approximated by a simple two-state model 



T 



Vb/16 

Vo/16 (q+1) 2 



(8) 



h q can be brought into the form of the Hamiltonian 
given by Eq. ([!} by properly shifting the diagonal parts 
(e.g., shifting away the quadratic term in time t ). For 
q = we thus immediately obtain: 



2F E Tec t /nh 

Vo/4 



Vb/4 
-2F Q E lec t /Trh 



(9) 



The a, AE and 7 introduced in the LZ model of 
Eq. (fTJ) can be expressed in terms of our system param- 
eters: a = 2F ElJjrh = 4E tcc /(ttT b ), AE = V E Tec /2, 
and 7 = 32Fq/ttV 2 . The LZ theory predicting the 



asymptotic behavior of the tunneling probability, can be 
used as a very good approximation for our system for 
times far enough from the avoided crossings. However 
there are some limiting cases, and experimental param- 
eters, for which the simplified two-state model is not 
a good approximation for the Wannier-Stark system. 
The discrepancy is large for lattice depths larger than 
the energy scale E lec of the system (Vo S> 1), where the 
gap between energy bands increases leading to quasi-flat 
bands and localized eigenstates. Therefore, several 
momentum eigenstates contribute with a non-negligible 
amount to the lowest energy eigenstate, and one would 
need to take into account more components in the 
Hamiltonian matrix. 




FIG. 4: Comparison between the time evolution of the Bose- 
Einstein condensate survival probability, in the (a) diabatic 
and (b) adiabatic basis, for different initial conditions, pre- 
pared at a temporal distance At = Tb/2 from the crossing 
point. The dashed line is the evolution in a possible experi- 
mental setup, i.e. the evolution following an initial prepara- 
tion in the ground state of the adiabatic basis; the solid line 
is the evolution following an initial preparation in the ground 
state of the diabatic basis. Parameters are Fq = 1.197 and 
Vo = 2.3, corresponding to 7 = 2.3, leading to a jump time in 
both adiabatic and diabatic bases 1.9 times the Bloch period 
T B . 



B. Initial conditions 

Before analyzing the experiment results, we need 
to address an additional problem, the finite coupling 
duration, as in [301 ] : an experiment necessarily takes a 
finite time for the measurement, whereas the standard 
LZ theory assumes that the time taken for the transition 
runs from —00 to 00. The experimental finiteness of the 
sweep time Tb implies that for the initial state at a finite 
distance from the transition point, the diagonal and 




FIG. 5: Comparison between the time evolution of the Bose- 
Einstein condensate adiabatic survival probability, starting 
from initial ground states prepared at different time distances 
from the transition point. Survival probabilities measured for 
F = 1.197 and V = 0.3 in (a), and for F = 1.197 and 
Vo = 3.0 in (b). The dotted lines show the evolution ob- 
tained evolving the survival probability from the ground state 
in the adiabatic basis at t — —00; the solid lines illustrate the 
evolution obtained evolving the survival probability from the 
ground state which simulates a possible experimental initial 
state, i.e., the ground state in an adiabatic basis, at a finite 
time from the transition point. 



off-diagonal matrix elements in the system Hamiltonian 
are comparable. The experiment we are dealing with 
typically operates in the regime defined in [30j | as a large 
time, meaning that the time intervals from the turn-on 
and the turn-off times to the crossing are larger than the 
jump time. The presence of a jump time comparable to 
the Bloch time may modify the temporal evolution of the 
survival probability for the two mechanisms discussed 
in the following. Because at large 7 the ratio T^ ump /Tb 
between jump time and sweep time is given by 7P\/7o/2, 
large Fq values may produce deviation from the ideal 
theory of [H. 

At t = —00 the diabatic and adiabatic states coincide, 
and hence, the preparation of the system in its ground 
state is unambiguous. On the other hand, at a finite 
distance from the transition point, the diabatic and adi- 
abatic states do not coincide. In a numerical approach 
any chosen initial state can be evolved given the proper 
Hamiltonian, but from an experimental point of view 
the system will be prepared in a well defined initial 
state, which depends both on the parameter values 
and on the preparation protocol. It is not obvious that 
this initial state can be chosen at will: most likely, the 
experimental initial state will be the one corresponding 
to the ground state of the complete Hamiltonian, i.e., 
the adiabatic lower energy state, at a time equal to the 



time when the sweep starts. The comparison between 
experiments and theory performed for different initial 
states should clarify this issue, because the evolution 
for different initial states is markedly different, both 
when observed in the diabatic basis and in the adiabatic 
basis, see Fig. EJa) and (o) corresponding to typical 
Bose-Einstein condensate experimental parameters. The 
results of Fig. @]show that, for experimentally accessible 
parameters, the two evolutions do not coincide in both 
the diabatic basis and adiabatic bases (see Fig. \M. a ) 
and (6), respectively). We have verified that the results 
of Fig. 3] following an initial preparation in the ground 
state of the diabatic basis (solid lines) coincide with the 
finite coupling duration predictions of ref. [30J ■ 

It is not at all obvious that an initial state chosen 
as the adiabatic ground state at a finite time from 
the transition point (which is likely to be the initial 
experimental state) should coincide with the state 
obtained evolving from t = — oo, projected onto the 
adiabatic basis. We computed the survival probability 
simulating different Bose-Einstein condensate initial 
states, see Fig. [5J For our experimental parameter set, 
the discrepancy is not very large, but certainly important 
for a precise description of the temporal evolution of the 
tunneling. Therefore, the approach of [29|] yields some 
elegant theoretical results for the LZ transition, but care 
is needed in comparing them with the experiment due 
to the presence of the additional time scale connected 
to the finite distance between the experimental starting 
point and the transition point. 



C. Numerical calculation 

In [381 ] some of us have introduced an easily computable 
quantity to determine in a good approximation the sur- 
vival probability in the adiabatic basis: 



/oo 



dp\^( P ,t)\' 



(10) 



where ^(p, t) is the Bose-Einstein condensate wave func- 
tion in momentum representation, and p c > 3p lcc is an 
ad hoc cut-off. Eq. (fTU|) can be interpreted as the projec- 
tion of ^(p, i) onto the support of the initially prepared 
condensate at t = (in the presence of the optical lattice 
but at Flz = 0), which is illustrated in Fig. [61(a). Since 
Eq. (jTU)) measures the decay only once the Bose-Einstein 
condensate wave packet ^(p, t) has extended beyond — p c 
(= — 'ip r cc in Fi g R)f a)), we must resort to the accelera- 
tion theorem |6j,|26[ to identify time t with £ — Tb, i.e., we 
must rescale time by the traversal time of the Brillouin 
zone Tb- 

While many previous experimental results proved in 
very good agreement with simulations based on Eq. (fTDj) . 
c.f. [34M36J. a better numerical method is needed for 
the new generation of experiments reported here. The 
dash-dotted lines in Fig. [5] were produced using P a (t) of 



Eq. ((TO)) . These simulations well reproduce the height of 
the steps in agreement with the LZ prediction given in 
Eq. ^. They do not, however, reproduce the oscillations 
of the experimentally measured survival probability, due 
to the artificial cut-off used for evaluating P a (t). While 
the sequence of steps - corresponding to a sequence of 
LZ tunneling events - is observable in Fig Et^»)i n o oscil- 
lations arc visible. To reproduce the oscillatory behavior 
of the experimental data in Fig. [2l instead of Eq. (fT0|) 
we determine Pn(t) in the following way: \(f>(n,q)) shall 
denote the band solution for the ground band n = as 
shown in the lower left panel of Fig [31 Then the adiabatic 
survival probability is just the projection of the conden- 
sate wave function ^(p, t) onto (f>(n = 0, q) integrated 
over the full Brillouin zone, i.e., 



P*(t) 



dq\(y(p = q,tM0,q))f 



(11) 



The survival probabilities P a (i) shown in figs. 2] and [5] 
have been calculated in this way. 
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FIG. 6: (Color online) In (a) the momentum distribution at 
time t = 10Tb (dashed lines) starting from the initial mo- 
mentum distribution (solid lines) under the action of a force 
directed towards negative p values. The vertical dash-dotted 
line shows the cut-off value p c — 3p ICC in the definition of 
Eq. (|10[) . (b) Temporal evolution of the survival probability 
in the adiabatic basis using mentioned definition. Simulation 
parameters: V = 2.07, F = 1.197. 

On the other hand, following the procedure sketched 
in the lower right panel of Fig [31 the survival probabil- 
ity determined in the diabatic basis of free momentum 
eigenstates is given by 



Pd(t) 



dq\(*(p = q,t)\p = q)\ 2 , (12) 



with p — q within the first Brillouin zone in the notation 
of Section IIA. Eq. ([T2"]) is used to simulate the experi- 



mental results of Figs. 0/0 and (b) presented in the next 
section. 



IV. RESULTS 

In our experiments we realized the Wannier-Stark 
Hamiltonian of Eq. ((U) with Bose-condensed rubidium 
atoms inside an optical lattice [13, |32h36| . Initially, 
we created Bose-Einstcin condensates of 5 x 10 4 87 Rb 
atoms inside an optical dipole trap (mean trap frequency 
around 80 Hz). A one-dimensional optical lattice created 
by two counter-propagating, linearly polarized gaussian 
beams was then superposed on the Bose-Einstein 
condensate by ramping up the power in the lattice 
beams in 100 ms. The wavelength of the lattice beams 
was A = 842 nm, leading to a sinusoidal potential 
with lattice constant d L = A/2 = 421 nm. A small 
frequency offset Ai/(i) between the two beams could 
be introduced through the acousto-optic modulators in 
the setup, which allowed us to accelerate the lattice in 
a controlled fashion and hence, in the rest-frame of the 
lattice, to subject the atoms to a force Flz = Maui 

.,, j dAv(t) 

with a LZ = d L dt y . 

In several previous experiments [32h35| we had 
already measured the LZ tunneling probability by first 
loading the Bose-Einstein condensate into a lattice, 
then accelerating the lattice for one Bloch period (i.e., 
across the zone edge and then to the center of the second 
Brillouin zone) and subsequently measuring the number 
of atoms left in the fundamental band. This was done 
by accelerating the lattice further with a smaller value 
of a scp and a larger lattice depth V scp chosen such as 
to ensure that atoms in the fundamental band did not 
undergo LZ tunneling at subsequent crossing of the zone 
edge and that atoms in higher bands tunneled with 
almost 100% probability. In that way it was possible to 
separate atoms in the fundamental band in momentum 
space so that after a time-of-flight they could be easily 
measured. 

The time-resolved measurements we are interested in 
for the purposes of the present paper initially followed 
the same procedure. Rather than accelerating the lattice 
for a full Bloch period, however, we had to interrupt the 
LZ tunneling event at some time t ^ nTe in general. The 
exact protocol then depended on whether we wanted to 
measure in the adiabatic or in the diabatic basis. 

For measurements in the adiabatic basis, we proceeded 
as follows, see Fig. [3J After loading the Bose-Einstein 
condensate into the optical lattice, the lattice was 
accelerated with acceleration <2lz for a time £lz- The 
lattice thus acquired a final velocity v = olz^lz- At 
time t = t-^z the acceleration was abruptly reduced to a 
smaller value a scp and the lattice depth was increased 
to V^ep in a time i ram p <C Tb- These values were chosen 
in such a way that at time t = £lz the probability for 
LZ tunneling from the lowest to the first excited energy 
band dropped from between at 0.1 — 0.9 (depending 



on the initial parameters chosen) to less than w 0.01, 
while the tunneling probability from the first excited to 
the second excited band remained high at about 0.95. 
This meant that at t = £lz the tunneling process was 
effectively interrupted and for t > i L z the measured 
survival probability P(t) = No/N to t (calculated from 
the number of atoms Nq in the lowest band and the 
total number of atoms in the condensate N iot ) reflected 
the instantaneous value P(t = £lz)- 

The lattice was then further accelerated for a time 
t sop such that a SO p£ sep ps 2mp lcc /M (where typically 
m = 2 or 3). In this way, atoms in the lowest band 
were accelerated to a final velocity v ps 2mp Tec /M, 
while atoms that had undergone tunneling to the first 
excited band before t = £lz underwent further tunneling 
to higher bands with a probability > 0.95 and were, 
therefore, no longer accelerated. At time t scp the lattice 
and dipole trap beams were suddenly switched off and 
the expanded atomic cloud was imaged after 23 ms. In 
these time-of-flight images the two velocity classes and 
2mp rcc /M were well separated and the atom numbers 
N and N tot could be measured directly. Since the 
populations were effectively "frozen" inside the energy 
bands of the lattice, which represent the adiabatic 
cigenstatcs of the total Hamiltonian of the system, 
this experiment measured the time dependence of the 
LZ survival probability P a in the adiabatic basis, c.f. 
Eq. (JTJJ) above. 

The results of our measurements in the adiabatic basis 
are summarized in Fig. [2] The step-like behavior of the 
survival probability around t = 0.5Tb is clearly visible, 
as well as the finite width of the step, which demon- 
strates that our experimental protocol does, indeed, 
allow us to access the dynamics of the LZ transition 
and the jump time associated with that transition. 
Also shown in the figure are the results of numerical 
simulations using the cut-off and the adiabatic survival 
methods described above in Section IIIC. As expected, 
both methods reproduce the step with a finite width 
and the steady-state value of the survival probability 
for long times. The slight oscillations of the survival 
probability for t > 0.5Tb, however, are only visible in 
the results computed according to Eq. (fTTj) above. In 
fact, the amplitude of these oscillations is larger in the 
numerical simulations than in our experimental data. 
This might indicate that our protocol for freezing the 
instantaneous populations in the ground and excited 
bands is not perfect. Indeed, we found that a delicate 
balance between the accelerations and lattice depths 
for the separation phase was necessary in order to 
ensure that the populations after the separation phase 
faithfully reproduced those at t = t^z, which was tested 
by choosing two extreme values for olz which gave 
theoretical survival probabilities of approximately and 
1, respectively, and then verifying that these values were 
measured in the experiment. In practice, the parameters 
for the separation phase were optimized in this way for 
one set of the LZ parameters and then kept constant as 



V was varied in Fig. [2j 

For measurements in the diabatic basis, the exper- 
imental protocol was even simpler, see Fig. |3] As in 
the adiabatic case, after the initial loading phase the 
lattice was accelerated with acceleration olz for a time 
tuz. At that point the atomic sample was projected 
onto the free-particle diabatic basis by instantaneously 
(within less than 1 fj,s) switching off the optical lattice. 
After a time-of-fiight the number of atoms in the v = 
and v = 2p rec /M momentum classes are measured and 
from these the survival probability (corresponding to the 
atoms remaining in the v — velocity class relative to 
the total atom number) is calculated. Fig. [7] shows the 
results of such measurements, together with numerical 
simulations based on Eq. (fl2|) . As in the adiabatic case, 
a step of the survival probability around t = 0.5 Tb is 
clearly seen, as well as strong oscillations for t > 0.5 Tb. 
These oscillations are much stronger and visible for a 
wider range of parameters in the diabatic basis than 
in the adiabatic basis (see the results for Vo = 2.3 in 
Fig. [2J which is confirmed by our numerical simulations). 
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FIG. 7: (Color online) Time-resolved measurements of LZ 
tunneling in the diabatic basis, a) Fixed force To = 1.197 
with different lattice depths Vo = 1 (filled circles), 1.8 (open 
circles), 2.3 (filled squares) and 4 (open squares), b) Fixed 
lattice depth Vo — 1.8 with different forces: To = 2.394 (open 
circles), 1.197 (filled squares) and 0.599 (open squares). The 
dashed lines are the results of numerical simulations based on 
Eq. ()12p which nicely reproduce the experimental data. 



Ultracold atoms in optical lattices provide an ideal 
model system for time-resolved studies of LZ tunneling. 
The complete control over the parameters of the lattice 
makes it possible to measure the tunneling dynamics in 
the adiabatic and diabatic bases by using different mea- 
surement methods. Our results confirm the existence of 
a finite temporal width for the transition in both bases 
and of strong oscillations of the survival probability in 
the diabatic basis. Both of these features arc backed up 
by numerical simulations taking into account details of 
the experimental protocol. 

Our findings pave the way towards more quantita- 
tive studies of the tunneling time for LZ transitions, 
which are of current interest in the context of optimal 
quantum control and the quantum speed limit [39|. 
Also, it should be possible to measure the tunneling 
dynamics in arbitrary bases by inducing a rotation of 
the 2x2 LZ-matrix through variations in the lattice 
depth during the transition. With an appropriate 
choice of this variation one could then, for instance, re- 
alize the super-adiabatic basis proposed by M. Berry [281 ] . 
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